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ABSTRACT 


We  analyze  the  convergence  of  the  spectral  vanishing  method  for  both  the 


spectral  and  pseudospectral  discretizations  of  the  inviscid  Burgers'  equa¬ 


tion.  We  prove  that  this  kind  of  vanishing  viscosity  is  responsible  for  a 


spectral  decay  of  those  Fourier  coefficients  located  toward  the  end  of  the 


computed  spectrum;  consequently,  the  discretization  error  is  shown  to  be  spec¬ 


trally  small  Independent  of  whether  the  underlying  solution  is  smooth  or 


not.  This  in  turn  implies  that  the  numerical  solution  remains  uniformly 


bounded  and  convergence  follows  by  compensated  compactness  arguments. 
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INTRODUCTION 


In  this  paper,  we  extend  the  analysis  of  the  spectral  vanishing  viscosity 
method  for  stabilizing  spectral  approximations  of  nonlinear  conservation 
laws.  The  spectral  vanishing  viscosity  has  been  first  introduced  in  [3], 

where  it  was  shown  that  L°° — bounded  spectral-Galerkin  approximations  converge 
2 

strongly  in  ^loc^^’*"^  exact  entropy  solutions  of  such  conservation 

laws . 

The  analysis  is  performed  on  the  2TT-periodic  inviscid  Burgers'  equation 


(1.1) 


-  h  ■ » 


submitted  to  the  additional  entropy  condition 


3_  <  0 


which  singles  out  the  unique  "physically  relevant"  weak  solution  of  (1.1). 
Both  the  spectral-Galerkin  and  pseudospectral-collocation  methods  for  (1.1), 
(1.2)  are  treated,  and  to  this  end  we  proceed  as  follows. 

Denote  by  S[,ju(x,t)  the  spectral-Fourier  projection  of  u(x,t). 


(1.3)  S  u(x,t)  = 


lk!^< 


u(k,t)e 


ikx 


u(k,t) 


2tt 


2Tr 


/ 

0 


u(x,t)e  ^'^^dx. 


and  lot  TJ^,u(x,t)  denote  the  pseudospectral-Fourler  projection  of  u(x,t), 
which  interpolate  u(x,t)  at  the  2N+1  equidistant  collocation  points 


IKX 
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(1.4)  I  u(x,t)  =  I  u(k,t)e  u(k,t)  = 


2N 


-ikx 


^  |k|  <  N  v=0  ^ 


These  two  projection  operators  differ  by  aliasing  error — that  is,  we  have 


(1.5) 


where  the  aliasing  projection  Am  is  given  by  [2] 


(1.6) 


A„u(x,t)  =  y  [  y  u(k  +  j(2N  +  1),  t)]e 
|k[<N  j^O 


ikx 


Throughout  this  paper,  we  use 


(1.7) 


as  a  concise  notation  for  the  two  kinds  of  Fourier  projections:  either  with 
a  =  0,  corresponding  to  the  spectral  projection,  or  with  a  =  1  which  corre¬ 
sponds  to  the  pseudospectral  interpolation. 

We  approximate  the  Fourier  projection  of  the  exact  solution  Pj^u(x,t),  by 
an  N-trigonometric  polynomial,  U|^(x,t) 

(1.8)  u„(x,t)  =  y  u,  (t)e^'^’^, 

^  !kT<N 

which  is  determined  by  the  approximate  evolution  equation 


(1.9)  u^(x.t)  +  (Pj^  i  u^  (x,t))  =  e 


3x 


3x 


(x,t)) 
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The  expression  on  the  right-hand  side  of  (1.9)  represents  the  spectral  vanish¬ 
ing  viscosity  term.  Here  is  the  spectral  viscosity  operator  which  is 
defined  as  a  convolution  with  a  symmetric  viscosity  kernel  Qj,j(x), 


k  <N 


In  the  spectral  case  where  a  =  0,  (1.9)  amounts  to 


It  It  T  =  e  [0^(x)*  Uj^(x.t)] 


J 


consisting  of  a  nonlinear  system  of  ordinary  differential  equations  for  the 
Fourier  coefficients,  Uj^(t),  which  are  coupled  through  the  standard  spectral 
convolution  treatment  of  the  nonlinear  term.  The  interpretation  of  the  scheme 
(1.9)  in  the  pseudospectral  case  where  a  =  1  leads  us  to 


It  It^^n  I  |x=x  = "  It  |x=x  .  o  <  v  <  2n, 

V  V 


and  consists  in  a  complete  statement  of  a  standard  collocation  method  with  a 
pseudospectral  treatment  of  the  nonlinear  term. 

In  both  the  spectral  and  pseudospectral  cases,  the  spectral  viscosity 
operator  can  be  efficiently  implemented  in  the  Fourier  rather  than  the  physi¬ 


cal  space ,  1 .e . , 


^  It  It  "  It  It  =  '"‘lki'<N 


An  essential  ingredient  of  our  spectral  viscosity  operator,  Q^j,  is  that  it 


I 


should  operate  only  on  the  high  portion  of  the  spectrum,  in  order  to  retain 


the  formal  spectral  accuracy  of  the  method.  Hence  we  make 


ASSUMPTION  I:  There  exists  a  constant 


m  =  m(N)  <  —  N, 


such  that 


(qCk)  E  0, 

'  0  <  Q(k)  <  1, 

I  Q(k)  =  1, 


|kl  <  m, 
m  _<  I  k  I  <  2m , 
2m  <  |k|  <  N. 


Then,  with  =  I  -  we  can  rewrite  (1.9)  as 


(1.11)  u^(x,t)  +  (P^  I  =  e  —  [(I  -  R^)  u^(x,t)]. 


where  the  corresponding  kernel  Rj^Cx), 


(1.12) 


R  (x)  =  y  R(k)e^'^^, 
|kT<2m 


Is  a  trigonometric  polynomial  of  degree  _<  2m,  with  Fourier  coefficients 

|k|  <  m, 
m  _<  |k  I  _<  2m. 

In  order  to  guarantee  the  uniform  boundedness  of  our  approximation,  u^(x,t). 


(1.13) 


(R(k)  i  1, 

^  0  <  R(k)  <  1, 


we  shall  need  to  control  the  size  of  this  kernel;  we  therefore  make 


-5- 


ASSUMPTION  II:  There  exists  a  constant  such  that 

(1.14)  IIR(»)ll,  <  Const.logm. 

m  1 ,  V  — 

L  (x) 

We  remark  that  the  assumption  of  a  logarithmic  upper  bound  for  the  size  of 

A 

R^(x)  is  plausible,  since  typical  applications  Involve  R(k)  which  de¬ 
crease  monotonically  to  zero  and  (1.14)  is  automatically  fulfilled  in  such 

cases;  consult  Appendix  A.  To  obtain,  with  the  help  of  Assumption  II,  the 

00 

promised  uniform  bound  on  Ufj(x,t),  necessitates  L  -bounded  initial  data, 
Ufj(x,0).  For  technical  reasons  we  shall  need  the  slightly  stronger 

ASSOMPTION  III;  There  exists  a  constant  such  that 

A 

llu,,(x,t=0)»  <  y  (u,  (t“0)|  <  Const^. 

N  L*(x)  -  lkt<N  ® 

The  spectral  viscosity  term  on  the  right  of  (1.11)  depends  on  two  free 
parameters:  the  viscosity  amplitude  e  =  e(N)  and  the  effective  size  of 

the  invlscld  spectrum  m  =  m(N).  These  two  parameters  should  be  chosen  to  en¬ 
sure  the  convergence  of  the  method.  In  [3]  it  has  been  proved  that  in  the  ab¬ 
sence  of  such  viscosity  terra,  e  =  0,  strong  as  well  as  weak  convergence  to 
the  exact  entropy  solution  falls. 

The  main  result  of  this  paper  asserts 

Theorea  1.1:  Consider  the  Fourier  approximation  (1.11)  of  either 

spectral  or  pseudospectral  type.  Let  the  spectral  viscosity  in  (1.12)  - 
(1.14)  be  parameterized  with  (e ,m)  as  follows 


(1.15)  e  = 


Then  Uj>j(x,t) 
conservation 


eCtn)  ~  a  — 75 -  ,  m  =  m(N)  ~  Const. N^,  0  <  6  <  -^  . 

m^^llR  (•)«  , 

“  h^x) 

converges  boundedly  a.e.  to  the  unique  entropy  solution  of  the 
law  (1.1). 


Let  us  examine  for  example  the  viscosity  operator  =  I  -  Here 

Rjj,(x)  coincides  with  Dirichlet  kernel  Djjj(x),  where  [5,  Chapter  II] 


m 


Ikx 


k  <m 


1 

2n 


sin(m  +  ■^)x 


sin  X 


IID  (  •  )  II  1 - K  logm, 

L  (x)  TT^ 


so  that  Assumption  II  is  fulfilled  and  Theorem  1,1  yields 


Corollary  1.2:  Consider  the  Fourier  approximation 

(1.16)  ^  Uj^(x,  t)  +  ^  (Pj^  Uj^(x,t))  “  e:  -^  [  (I  -  S^)  Ujj(x,t)]  , 


wl  th 


(1.17)  e 


e(N)  ~  Const. 


logN  ’ 


m  =  m(N)  ~  Const. N^,  0  <  3  ^ 


Then  unj(x,t)  converges  boundedly  a.e.  to  the  unique  entropy  solution  of  the 
conservation  law  (1.1) . 


The  spectral  portion  of  this  result  (a  =0),  was  derived  in  [3,  Theorem 
4.1]  under  the  assumption  that  the  numerical  solution  U]>j(x,t)  remains  uni¬ 
formly  bounded.  The  extension  of  Corollary  1.2  includes  the  pseudospectral 


approximation  (a  =  1),  and  in  addition,  thanks  to  the  slightly  more  stringent 
parametrization  than  that  of  [3,  Theorem  4.1],  contains  a  proof  of  the  previ¬ 


ously  assumed  L  -bound. 

In  the  last  example  the  viscosity  symbols  Q(k)  were  discontinuous 
at  |k|  =  m.  It  was  suggested  in  [3]  that  the  use  of  viscosity  operators 
0,,  with  smoothly  varying  symbols  would  be  advantageous  for  the  spectral  vis¬ 
cosity  method  in  (1.9).  As  our  second  and  final  example  we  consider  the 
simplest  viscosity  operator  of  this  type,  namely 


|k|  <  m, 
m  ^  Ik|  _<  2m, 
2m  <  |k|  <  N. 


This  kind  of  spectral  viscosity  is  intimately  related  to  the  Fejer  operator 
j  m-1 

F  =  _  •  ^  S,  :  if  we  let  ^m(x)  denote  the  corresponding  Fejer  kernel  [3, 

k=0 

Chapter  III] 


K 

m 


(x) 


then  for  =  I  - 

associated  with 


(1  - 

|kl .  ikx 

=  ^  ( 

•  ^  9 

Sin  mx  2 

7  1 

IIK  (Oil  , 

”  l'(k) 

■  ■  •  /  ~ 
m 

TTin  ^ 

I  ^  ’ 

sin  Y  X 

^m 

we  have 

Rm(x) 

=  2K2„(x)  - 

K^(x).  Hence 

R 

m 


2F 


2m 


,  2m- 1 

F  =  -  •  I 

m  m  ,  ^  k 
k=m 


=  1, 


the  kernel 


i  5 


-uni formly 


bounded  , 


HR  (  •  )  II 
m 


(x) 


2IIK-  (.)ll  , 

2"^  L*(x) 


+  II K  (Oil 

LOx) 


<  3, 
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{?V^ 


so  that  Assumption  II  is  fulfilled  and  Theorem  1,1  yields 


Corollary  1,3:  Consider  the  Fourier  approximation 


«  ^111 
[(I  •  I 

T«  t-' 


2m- 1 


k=m 


with 


(1.19)  e  =  e(N)  ~  Const, N  ,  m  =  m(N)  ~  Const, N^,  0  <  g  <  -I 


Then  ujq(x,t)  converges  boundedly  a.e,  to  the  unique  entropy  solution  of  the 
conservation  law  (1,1), 


The  paper  is  organized  as  follows.  In  Section  2  we  derive  a  couple  of 
basic  L  -type  a'priorl  energy  estimates.  In  Section  3,  these  estimates  are 
used  in  order  to  study  the  spectral  decay  rate  of  the  Fourier  coefficients. 
This  in  turn  enables  us,  in  Section  4,  to  obtain  l“  a  priori  estimate  on 
the  numerical  solution.  Finally,  based  on  the  a  priori  estimates  prepared  in 
Sections  2,  3,  and  4,  Theorem  I.l  is  proved  in  Section  5  along  the  lines  of 
[3],  using  compensated  compactness  arguments. 


2.  L^-TYPE  A  PRIORI  ESTIMATES 

We  consider  the  approximate  Fourier  method  (1.9)  which  we  rewrite  as 


o  1  ^  9  9  /I  2,,  3 

^  ^  3t  “n  3x  ^2 


In  order  to  prove  the  convergence  of  this  method  we  need  a  couple  of  a  priori 
estimates  on  its  solution.  To  this  end,  we  multiply  (2.1)  by  Uj^j 


^  ,1  2.  3  A  3,  _ 

3t  2  3x  ^3 


(2.2) 


3x  2  '''  3x  3x 


and  integrate  over  a  2-n  period:  the  integral  of  the  second  term  on  the  left 
vanishes  by  periodicity,  and  after  integration  by  parts  for  the  second  term  on 
the  right  we  are  left  with 


(2.3) 


I  +  e  / 

L  (x)  0 


27T 


u^(x,t)Q^  u^(x,t)dx 


/ 

0 


2tt 


Using  (1.7)  and  the  fact  that  I  -  is  orthogonal  to  our  N-space,  we  find 

that  the  right-hand  side  of  (2.3)  equals 


2n 


"SIT  f 


1 

7 


Ztt 


u.Jdx  =  -a- 
N 


/ 

0 


P|<N 


^ip(A^  4 


N'p 


and  by  the  aliasing  relation  (1.6),  this  does  not  exceed 


-5^[(i-Pn)4 


u^ldx 


3  X  .  "  "  "  .  |a|N  7  ^ 

nr  ipU  U  U  < 

|p+q+rf=2N+l  P  ^  ^  |p+q+r|=2N+l 


u 

P 


In  view  of  |p  +  q  +  r 
! q I  <  N  and  I r |  <  N 


=  2N  +  1,  at  least  two  of  the  three  indices 

M 

are  greater  in  absolute  value  than 


|pi  <  N’, 


N 

1  ' 


and  hence 


J  “n  1  nr"  !pMu  iiq 

”  |<|p|<N  |<|q|<N 


'^q'  '^±(2N+l)-(p+q) 


E 

^  I  p I <N 


9  .  ^  ,9 

p“|u 


|2  2^  |2- 

p'  „2^  "  - 

|<|q!<S 


•[  E  E 

XT  XT 


^|p|<N  -^<|q|<N 


±(2N+1 )-(p+q) 


Consequently,  since  Q(p)  =  Q(q)  =  1  for  jp|,|q|  >.  T  ’  expression  on 


the  right  of  (2.3)  can  be  upper  bounded  by 


0  N  L  (x)  L  (x) 


'loreover,  since  0  <  0(k.)  <1  we  have  for  the  second  term  on  the  left  of 


(2.3) 


(2.4b) 


]  -5—  u,.(x,t)0,,  -5—  u.,(x,t)dx 

U  dxN  N3xN 


=  e  \  Q(k)|u^(t)|^  >  cIlOj^  Uj^(*  ,t)ll^, 

!k!<N’  ■  L‘'(x) 


Inserting  this  together  with  (2.4a)  into  (2.3)  we  end  up  with 


(2.r  T^"r('n)ii%  +  [e  -  l[^iiu  (•  ,t)!i  ,  ].iio  ■^u^^^(-,t)ii  ,  <  o. 

“  ■'  L^x)  ^  r(x)  ^  ^  L^x) 


Thus ,  ns  1 ong  ns 


A  q  I  c 

--fTn  llo  (  •  ,t)  il  ^  >  T  ' 

r}'~  T/ixi 


1  h  r  n  i  n 


In  particular,  (2.7)  implies  that  for 


we  have 


u»(x,t)  *  I  u 

,  IM<n  " 

our  first  L^-type  a'priori  estimate 


(2.8)  .  5  <  E^,  <  Co„st„. 

L  (x)  |k|_<N  L  (x) 

Hence  (2.6),  (2.7)  and  consequently  (2.8)  prevail  for  all  time  provided  (2.6) 
is  valid  at  t  =  0,  i.e.,  we  require  that  in  the  pseudospectral  case  where  a 
=  1  V7e  shall  have 


(2.9) 


e(N)  > 


indeed.  Assumption  II  tells  us  that  this  requirement  is  fulfilled,  at  least 
for  sufficiently  large  N,  for 


(2.10) 


e  >  Const . 


^  >  8.Er,*N 

logN  0 


-1/2 


2B  <  y 


Furthermore,  temporal  integration  of  (2.7)  then  gives  us  the  second  a'priori 
estimate 


(2.11)  elOjj 

hoc**’'* 


=  e  /  |Q(k)u  (t)  I  “^dt  < 

t  m<  I  k  I  <N 


Jq  <  Constp. 


;rtM 


The  first  integral  on  the  right  does  not  exceed 


(3.4) 


<  2  11(1  S2^)gj^  Uj^(«,t)il  2  *11(1  2.  ' 

L  (x)  L  (x) 


In  order  to  estimate  the  second  term  of  the  last  product,  we  will  make  use  of 


the  following  lemma  whose  proof  is  postponed  to  the  end  of  this  section. 


Lemma  3.1:  Let  f.,  =  f.,(x)  and  a  =  g.,(x)  be  a  couple  of  N- 

-  N  N  -  °N  “N  - *■ - 


trigonometric  polynomials.  Then,  for  any  0  <  2k  <  N  we  have 


(3.5) 


1—  ["fN''  2.  %".2.  /  2,  ^n''-2.  J 

/k  L  (x)  L  (x)  L  (x)  L  (x) 


Lemma  3.1  with  f.,(’)  =  g„(*)  =  u,(»,t)  implies 
N  N  N 


(3.6)11(1  S2^)PjjU^(*,t)ll  2  1  2.  ^k^'Sx  ”  2.  v' 

L  (x)  /k  L  (x)  L  (x) 


Equipped  with  (3.4),  (3.6),  and  (2.8)  we  return  to  (3.3)  to  find  that 


(3.7) 


7  dt  ^2k^^N''’^^"  2.  .  i  ^2k^3x  2,  / 

L  (x)  L  (x) 


"  ^  "o"«  -  52k>  llT  2,  •'«  -  2,  , 

✓k  L  (x)  L  (x) 


i 
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i 


r>.v,' 


which  brings  us  to 


Theorem  3.2j  For  any  integer  s  >  0  there  exists  a  constant  C 


ConstCsjEg),  such  that  for  sufficiently  large  N,  N  >  2®»4in,  we  have 

1  ^s 


(3.8)  11(1  -  S^)u^(..t)ll  <  C  +  (1 

Y  L  (x)  e/N  c/n 


Proof :  Let  E|^(t)  abbreviates  the  quantity 


1  'i-(s+l)  -4®*eN^t 


]. 


(3.9) 


E|^(t)  -  11(1  2 

L  (x) 


In  view  of  the  inverse  inequalities 


2kE2k(t)  ^  11(1  -  S^k)  '^N^* 2.  .1  NE2k(t), 

L  (x) 


it  follows  from  (3.7)  that  Ej^(t)  satisfy 


(3.10) 


H  9  2E„*N^ 

—  E,,  (t)  <  -4ek^E„.  (t)  +  -  Ek(t), 


dt  2k  - 


2k 


m  <  k  <  N. 


/k 


Temporal  integration  yields  that  for  any  0  <  tQ  ^  t  we  have 


(3.1.)  l"  .  e 


-4ek^( t-t  - ) 


/k  T=t 


0 


and  therefore 


E„-B^  -4Ek^(t-t  ) 

'=2k")  <  --j^r-T  •  "f  "  •'=2k<'o' 

2e/k*k 


(3.12) 


The  a'priorl  estimate  (2.8)  implies  that 


max  E  (t)  <  max  E  (x)  <  E 
0<T<t  0<T<t 


and  in  view  of  (3.12)  we  have 


(3.13) 


2k 


E  .N 

<7^ 


2e/k*  k 


2 

-4ek  t  - 

Eo  ®  *Eo> 


If  we  choose  ~  \  (3.12)  we  find  that 


Eq-N" 

(3.14)  E^j^(t)  £  - ^ ^  •  max 


2e/k»k  t,  .. 


E^(t)  +  e  ^  ^ 


and  following  Krelss  [1],  we  can  use  this  to  improve  our  estimate  (3.13). 
Namely,  for  k  >  2m  we  can  use  (3.13)^  *=0  upper  bound  the  terms  max  E,^(t) 

T 

and  ^2k^'7^  right  of  (3.14),  and  obtain  the  improved  bound 


8E„*N^  2  8E„*N^  _  2 

-Eo-  >^>2- 

e/k*k  e/k*k 


Now  we  can  repeat  this  process,  and  by  induction  we  obtain  that  for  k  >  2®»m 
we  have 


8®E„*N^  s+1 


(3.15) 


E2k(t)  £ 


0 


8®E„*N^  s  ,-s+l 


e/k*  k 


-)  -En  +  (1  + 


^0 


e/T^*k 


-)  e 


-4 


•ek  t 


^0‘ 


Verification  of  the  induction  step  is  left  to  the  Appendix.  Finally,  (3.15) 

N  s 

implies  that  for  sufficiently  large  k  =  -r-  >  2  ^m,  we  have 


11(1  -  S^)u^(.,t)ll  2  =  < 


L^(x)  j 


32*8  .E.  s+1  32.8  •£„  s  ,-(s+l)  „2 

y  f  O'!  n  X  z',  .  0-,  -4  "  .eN  t  ^ 

<  [ - - - J  -E  +  (1  +  - J  e  .E 

e/N  ^  e/N  ^ 


and  3.8  follows. 

Our  parametrizaclon  in  ( I . 14) , ( 1 . 1 5)  implies  that  for  sufficiently 
large  N  we  have 


,-(s+l)  .2  ^  ^  N 

4  .eN  t  >  Const.  - r- 


2-2B 


, s+1 I T 

4  .logN 


't  >  N^^^.t, 


0  <  28  < 


as  well  as 


1  logN  1  1 

— —  <_  Const.  — I -  <  —  ,  0  <  Y  -y  “  26  , 

e/N  (-i  -26)  “ 

N 

and  Theorem  3.2  tells  us  that 

Corollary  3»Sl  For  any  integer  s  >  0  there  exists  a  constant  Cg 
such  that 


3/2 

(3.16)  11(1  -  Sj^)Uj^(.  ,t)ll  2  1  Cg*(N"®  +  e"^ 

-J-  L  (x) 

Corollary  3.3  indicates  the  spectral  decay  of  the  Fourier  coefficients 
|u(k)|  with  wavenumbers  |k|  2.  Y  which  in  turn  implies  a  similar  decay 


for  the  discretization  error,  I,  on  the  right  of  (3.1). 


For  the  latter  we 
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(3.17)  ll(I-Pj^)i  u^,(.,t)ll  2  =  ll(I-Sj^)^  u^(.,t)ll^2  +  a^llA^  ^  uj(.,t)ll^2 

L  (x)  L^(x)  L“(x) 


The  Fourier  coefficients  of  the  two  expressions  on  the  right  are  given 


respectively  by 


^  |k|  >  N, 

p+q-k=0  ^  ^ 


In  both  cases,  either  |p|  >  or  |q|  >  ;  hence  each  one  of  these 


coefficients  can  be  upper  bounded  in  a  standard  fashion  to  yield 


(I  ,t)ll  2  +a  II u^(«  ,t)ll  2  < 

L(x)  L(x) 


I  (l+a‘')*[ 


N< I kl<2N 


^  |u  (t)!^]  <  4E2.N.1I(I-S^j)u^(.,t)l|2 


|p|>7 


L'Cx) 


and  bv  (3.17)  this  is  the  same  as 


(3.18) 


<  2F,„./S  ..(I  -  »„)»/•. tX  ,  . 

L  (x)  ^  L  (x) 


Corollary  3.3  together  with  3.18  show  that  due  to  the  presence  of  the  spectral 


viscosity  term  IT  on  the  right  of  (2.1),  the  discretization  error  I  decays  to 


zero  at  a  spectral  rate  Independently  whether  the  underlying  solution  is 


smooth  or  not.  We  state  this  as 


Corollary  3.4:  For  any  integer  s  >  0  there  exists  a  constant 
such  that  for  sufficiently  large  N  we  have 


_  3/2 

(3.19)  III  :  [(i-P  )i  u;(*,t))»  <  C  ./N  -[n"®  +  e“^  s  >  n. 

^  ^  L^(x)  ~  ® 


V.'e  close  this  section  with  the  promised 


Proof  (of  Lemma  3.1):  Starting  with  the  identity 


and  subtracting  from  this  (I  -  H  (I  -  f|^  •  Sj^g^ )  =  0, 

we  can  write 


The  quantity  inside  the  right  brackets  is  a  trigonometric  polynomial  of  degree 
<  2N  and  hence,  by  Parseval  relation,  its  L^(x)  norm  dominates  the  l2(x) 
norm  of  its  Pv-  projection,  i.e.. 


(3.20) 


L  (x)  L  (x) 


1  ,• 


L  (x) 


The  norm  on  the  right  of  (3.20)  Is  upper  bounded  by 
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L  (x) 

(1.21  ) 

Ilf  II  ,  .ll(I-S  )£;  II  +  llg  II  .Il(l-S  )f  II 

L“(x)  L  (x)  L“(x)  L  (x) 


Finally,  for  hy,  equals  either  fj,]  or  g],j  we  have 


II  (I-S,  )h,.ll  < 

k  N  ■■ 


»  i  - 

I.  (x)  Ip|>k  P 


(1.22) 


|p|>k  p  IPI^k  ^  /k  L'^(x) 


and  (3.1)  follows  from  (3.20),  (3.21),  and  (3.22). 


4.  L  A  PRIORI  ESTIMATE 

The  classical  energy  method  can  be  used  to  show  that  the  solution  of 

(2.1)  remains  uniformly  bounded  during  a  small  finite  time.  The  method 

reflects  the  fact  that  for  sufficiently  smooth  initial  data,  say  with 
3‘- 


^  u  (x,t=n)  which  are  L^-bounded,  the  process  of  shock  formation  takes  a 
3x' 

finite  time,  during  which  -5 —  u,,(x,t)  remains  uniformly  bounded  and  a 

o  X  N 

couple  of  Sobolev  norms  could  be  a'priori  estimated  during  that  time. 

For  a  brief  initial  time  intervals,  we  can  do  better  with  regard  to  the 
smoothness  of  the  Initial  data,  as  told  by 


Lemma  4.1:  Consider  the  Fourier  approximation  (1.11)  -  (I.IA)  with 

Initial  data  Uj^(x,t  =  n)  such  that  Assumption  III  holds,  i.e.. 


(4. 1 ) 


y  Ui,  (  t=0  )  I  <  Const  „  . 
lkT<N  "  ^  " 


Then  for  t  <  --  we  h av e 


(4.2) 


llu.^,(  •  ,  t )  II 


<  2 


C  (  x) 


I  kT<N 


U|^(t=n)| , 


t  < 


8Constp*N 


Proof:  Tne  Fourier  transform  of  (1.9)  reads 


(4.3) 

u  (t)  +  ik[  \  u  (t)u  (t)  +  a»  *  ^  u  (t)u  (t)]  =  -ek^Q(  k)u  ( t) . 

Multiply  the  real  (imaginary)  part  of  this  by  sgn(Reu^(t))  (respectively, 
s,gn(ImUj^(  t) ) )  and  sum  over  all  k's:  since  the  right-hand  side  is  negative  we 
obtain  after  summing  both  parts 


d 

dt 


I  k[<N 


|u  (t)|  <  (1  +  |a|)*2N-  y  |u  (t)|*  y  |u  (t)|  < 

"  (pT<n  P  k  p 


(4.4) 


<  4N*( 


I  'V'" 

:1<M  ^ 


)|)^ 


Integration  In  time  yields 


00  1  I  |u.(t)l  <  - 

L  (x)  |kI<N  l-4Nt- 


k  <N 


|U|^(t=0)|  lkT<N 


u,  ( t=0)  I  , 


and  (4.2)  foil ows . 


To  obtain  L  bound  for  later  time,  we  shall  carefully  iterate  on 

the  LP(x)  norms  of  Uj^(x,t).  To  this  end,  we  multiply  (2.1)  by  pu^  ^ 

and  integrate  over  the  2Tr -period,  obtaining 


|P  + 


^  lP(x)  P^^  "" 


(4.5) 


p'l  p*/  k 


By  Corollary  3.4,  the  discretization  error  is  negligibly  small:  using  (3.19) 

1  2 

and  the  fact  that  (I  -  P.,)  u.,  is  a  trigonometric  polynomial  of  degree 

N  2  N 

<  2N,  we  have  for  any  s  >  0 


(4.6)  ll|_  ((I-Pf^4  p.  .  i 

L'^(x) 

3/2 

/m  .11^  [(I-Pj^)|  U^(.,t)]ll  2  <  CgN^.(N"%e"^  *S. 

L  (x) 


Therefore,  by  Holder  inequality,  the  first  integral  on  the  right  of  (4.5)  does 
not  exceed 


2  7T 

(4.7) 


Lq(x)  ^  2  ^  ^  lP(x)- 


<  p*  llu.,(*  ,t)  l|P  '  •€  .(N^”®+N^e  ^  '’^),  1  +  i  =  1. 

-  ^  lP(x)  ^  P  ^ 


The  second  integral  on  the  right  of  (4.5),  with  Ojjj  =  I  -  R^,  equals 


P*/  =  -ep(p-l)-/  u5""(^)  dx  -  tp./  uj  ^  °  (R^u^)dx 


2x  „  8u 
r  ..P“2,  N. 


2tt  ,  -  2 
■  p- 1  9 


r  r  ^  w  ■  w^'  wvw  vwv^ip  wn^iii  i»^iwji  im^yiMnw'^wr 


9- 


? 

^ 


The  first  term  on  the  right  hnnd  side  is  negative  for  any  even  Integer 
p  >  2;  for  the  second  term  we  use  Holder  inequality  as  before,  obtaining 


'  IT  Z  Z 

(4.8a)  ep./  u^  ^  L-  (R^u^)dx  <  e  pi!  u^^,  (  •  ,  t )  !)  ^  ^  *11— 


n 


3  x“ 


(x)  3x" 


L^(x) 


Now,  since  R  u.,  =  R  (x)*u.,(x,t)  is  a  trigonometric  polynomial  of  degree 
m  N  m  N 

<  2m,  (consult  (1.12)),  we  can  estimate  the  lP(x)  norm  of  its  derivatives 
as  follows  [5,  Chapter  X] 


(4.8b) 


il^(R  u„)ll  <  (2m)^*IIR  ,t)ll  < 

3x“  ”  ""  lP(x)-  ""  lP(x)' 


<  4m  HR  (On  ,  •  llu„,(‘  ,t)n 

L^x)  ^  lP(x) 


Using  (4.8a)  and  (4.8b)  we  conclude  that  the  second  integral  on  the  right  of 
(4.5)  is  upper  bounded  by 


(4.9)  p*/  u^,  f'^M  I""  ^  4p»em^«IIR  (Oil  ,  Mlu  (*,1)11^ 

J,  N  3x  N  3x  N  -  m  ^^(x)  ^  lP(x) 


Ve  recall  that  according  to  our  parametrization  (1.15),  ( • )  II  ,  <  a. 

Lhx)“ 

Hence,  equipped  with  (4.7)  and  (4.9)  we  return  to  (4.5)  to  find  that 


pllu  (•  ,t)  ll'’  ^  ••^llu^  (•  ,t)ll  < 

"  l’'u)  ^  iTcx)- 


,8/2 


p  (I  u  ( •  ,  t )  11  ^  ^  [  4a»  llu  (  •  ,  t )  II  +  C  ‘(N^  ^  +  N‘'e  *”  )  ]  ' 

i\  .Dy\  »P/\  S 


L^x) 


lP(x) 


or,  after  division  bv  the  common  factor  of  p  II  u  ( •  ,  t  )  II ^  , 

hP(x) 


id*  .  d  t  I  f  111  I  ~t  d  ■  I  ii' t  I*  '  ~  r  It  '  hit'-t’j'- 


~T~  i!  t  )  II  <  4a«  llii  r*  ,t)  II  +  C  -(N*  ‘  +N“e  '  ). 

LP(xr  iP(x)  " 

Finally,  we  integrate  in  time  obtaining  by  Cronwall's  inequality  that  fi-ir 
any  0  £  t^  _<  t , 


(4.10) 

4a(t-t  ) 

llu  (•,t)ll  <  e  •[!!□  (".t  )ll  +  C  ‘(N'"  ^•(t-t„)+  /N  • 

L'’(xr  ^  "  L^Cx)  *  " 

Letting  p  even  tends  to  infinity,  then  (4.10)  with  t^  =  *-0^'’^ 
gives  us 


e 


i/2 


•  t 


1 

RCons  1 *  N 


A  ^  t  r 

(4.11)  llu  (.,t)ll  ^  <e  .[llu^^(.,t  )ll  ^ 

L  (x)  L  (x) 


+  C  (N'"".t  +  /X 


and  together  with  Lemma  4.1  we  conclude  with  the  desired  L  bound,  namely, 


Theorem  4.2:  Consider  the  Fourier  approximation  ( 1 .  1  I )- ( 1  .  1 4  ) .  Then  for 

anv  s  >  0  there  exist  constants  a  >  0  and  C  such  that 
— ^  —  -  -  s - 


(4.12) 

^ CD  *'  e^^^*[4»  y  lu^(t=0)|  +  C  .^’(N  ^'t  +  /TT  e  j  ^ 

l'”(x)~  !kl£N  ^  ^  “ 

Remarks:  1.  We  observe  that  the  exponential  time  growth  in  (4.12)  does 

~  em“*ltR  f*)!!  , 


not  exceed 


4a  where  a 


<  Const. 


Kl-.- 


.  The  I'priori  Lf’(x)  estimate  deriv'od  in  (4, 10)  is  valid  for  .irbitrarv  L-- 

OO 

h'cod  initial  data.  We  note,  however,  that  the  resultinit  L  bound  in  such 
case  is  not  imiforin  with  respect  to  the  initial  time  tp.  That  is,  with  arbi- 

may  still  grow 


trarv  L“-bounded  initial  data,  the  solution  II  Uj^.l  •  ,  t )  II 


L  (x ) 

like  The  point  made  in  Lemma  4,1  was  that  with  slightly  strengthened 


regularity  assumption  on  the  initial  data 


kf<N 


I  u  ( t=0  )  I  <  Const  , 


0' 


t  hi 

on  . 

this  growth  is  bounded  for  a  brief  time  interval  of  length 


1 

N  ’ 


after  which 


5.  CONVERGENCE  TO  THE  ENTROPY  SOLUTION 

Wc  follow  [3],  using  compensated  compactness  arguments  to  conclude  that 
converges  to  the  entropy  solution  of  (1.1),  (1.2). 


Proof:  (of  Theorem  1.1).  Consider  the  four  terms  on  the  right-hand  side 


of  (2.1)  and  (2.2).  Using  (3.18)  together  with  (2.11)  along  the  lines  of  [3, 
Lemma  3.1],  we  find  that  the  term  I  satisfies 


(S.  !) 


< 


1 


<  ll(I-P^)^u^(-,t)ll  „  <  4F  .J  • 

L?  (X,t)  -  "  ^  N- 

1  oc 


->•  0. 


Also,  bv  the  a'priori  estimate  (2.11)  we  have 


*  '"'n  "n"  ?  ,  ,  '’• 

1 . ,  (  X  ,  t  )  N 

1  (^(' 


I’hus,  ill  view  o(  (^1.1)  ,111(1  ( 1 .  i’ )  ,  the  terms  f  .liul  II  on  t  lie  rlitlit  of  (2.1)  be 


1  oiiit  lo  the  I'omn.iel  of  II,  (x,l) 


Next  we  noti'  tli.it  sIikm'  (1^^  t-  I  we  li.ivi' 


'k  "n",  1  ,  <  "’n  r.  "n",  1  ,  "K  3  ,  , 


The  first  t »' rm  on  the  rl>.tht  is  bounded  by  .1^;  the  seeotid  one — be  i  rii'  tiie 


d  (' r  i  v.'i  r  i  ve  of  .1  t  r  i  ftotiomi't  r  i  r  polynomial  of  de^tfi'i'  <  2ni,  does  not  excei'd 


/i:  •  2m  ‘llti  II  ^  <  Const.  Consequently, 

I.y  (x,t) 

1  nr 


(  b  . '! ) 


^  ll’T—  II  o  Const  . 

5x  N  2  , 

h,  (x,t) 

1  or 


equipped  with  Cl.!)  we  now  turn  to  consider  the  right-hand  side  of  (2.2).  For 


the  third  term  in  (2.2), 


"n  ^  i.  "n' 


('■1.4a) 


f  u  ( I  -  P  )  i  u^l 
J  X  N  N  2  N 


SV  "n 


''n>  T  "n  ^  "'l  ^  "C’ 


we  hav('  by  Theorem  4,2  and  the  estimate  C).!) 


iSKcvKSiSJSc'SKKi 


( s,;b) 


II  u,.  (I  -  i\.ll  <  II  u^. 

I."  (x.t) 

1  or 


and  together  with  (5.3)  we  also  have 


<  llu  (•,t)ll  •11(1  -  P.,)-^  u,.ll  , 


I-  (x) 


I  “  (  X  ,  t  ) 

1  or 


mil.,  i  a^.  •  (I  -  P^,)^  u^ll 


(5.4c) 


<  /e  11^  U.,ll  ^  •  — ll(l-p  >4  T  <  flonst.4F  J  •  - ^ 

loc  loc 


Finally,  for  the  fourth  term  in  (2.2), 


(5.5a) 


57  ^^.N’  57 


3  r  A  9  ) 

^  57  57  “n’ 


^  A  ^ 

57  57  “n 


IV,  +  ns, 


have  by  (2.8),  (2.11)  and  the  uniform  bound  in  Theorem  4.2, 


(5Ab)  .rv,  c  S  (“s.",,  S  “n"„-1  ^,1  ■  '• 

loc’ 


while  TV 


2  =-  -=IS  S  "s-”s  S  "v 


CiAc)  IIV  .  <AO,,SvT  ,  , 

hoc'’'’'’  '■loc'’'-'* 


*  37  ? 


'447 


N’  57^'^2  . 


<  J“  +  /e  •2m*Ep«jQ  <  Const. 


Therefore,  by  Murat's  lemma,  [4],  the  inequalities  (5.4),  (5.5)  imply  that  the 

terms  111  and  IV  are  also  in  the  compact  of  H^^^(x,t).  In  summary,  we  have 

showTi  that  the  right  hand  sides  of  (2.1),  (2.2)  lie  in  the  compact  of 

H,  ^  (x,t),  and,  according  to  Theorem  4.2,  that  llu„(*,t)ll  is  bounded 

loc 

b  (x) 

(in  fact,  llu,,ll  with  p  >  6  will  do  for  our  purpose).  Hence  we  can 

”  Lhx.t) 

apply  the  div-curl  lemma  [4]  to  the  left-hand  sides  of  (2.1),  (2.2)  and  obtain 

9 

that  (a  subsequence  of)  UT,[(x,t)  converges  strongly  in  ^  weak 

limit  solution  u(x,t). 

Moreover,  we  claim  that  this  limit  is  the  entropy  solution  of  (1.1).  To 
verify  this  claim  we  show  that  the  right-hand  side  of  (2.2),  III  +  IV,  tends 
weakly  to  a  negative  measure.  Indeed,  by  (5.4)  and  (5.5b)  the  terms  III  and 
IVj  tend  weakly  to  zero,  and  hence  it  is  left  to  show  that  the  term  IV2> 


^  ^m  = 


tends  weakly  to  a  negative  measure.  To  this  end  we  proceed  as  in  [3,  Section 
4]  and  rewrite  IV2  in  the  form 


(5.6)  TV. 


rn  9  l2 

"  ^ 


+  e 


9x 


(R  u  )■ 
m  N 


Denote  the  three  terms  on  the  right  of  (5.6)  by  1^21.  ^^22  ^^^23’ 

respectively.  by  (2.11),  IV-,,  tends  weakly  to  a  negative  measure 


wlim  [IV-, 


a 


N  ^ 


<  0. 


(5.7-9) 


=  -e[Q 


-29- 


From  (5.7a)  -  (5.7d)  we  conclude  that  the  term  IV2  In  (5.6) — and  therefore 
that  the  right-hand  side  of  (2.2),  tends  weakly  to  a  negative  measure.  Thus, 

by  taking  the  weak  limit  of  (2.2)  we  recover  (1.2)  for  our  limit  solution 

_  2  — 

u(x,t).  Consequently,  the  strong  limit  of  Uj^(x,t)  =  u(x,t)  is 

the  unique  entropy  solution  of  (1.1)  as  asserted. 
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APPENDIX 


A.  THF.  L^-LOGARITHMIC  BOUND  OF  MONOTONF  VISCOSITY  KFRNFLS 

We  consider  symmetric  viscosity  kernels  of  the  form 

Qj^(x)  =  y  Q(k)e  +  y  e  , 

I  k|_<2m  2m<Tk  |  <N 

with  monotonical ly  increasing  Fourier  coefficients.  Then,  the  kernels  which 
correspond  to  Rj^  =  I  -  are  symmetric  polynomials  of  degree  <  2m 

(a.l)  R  (x)  =2*  y  R(k)coskx 

"  |kf<2m 

whose  Fourier  coefficients  are  monotonically  decreasing,  compare  (1.13), 

(a. 2)  1  2  R(k)  +20. 

Such  kernels  satisfy  Assumption  II  above,  as  told  bv 

Lenna  A.  1 :  There  exists  a  constant  such  that 

(a.l)  IIP  (•)I1  <  Const.logm. 

TT)  7  /  \ 

L  (  X  ) 

Proof;  The  result  follows  if  we  can  show  that  R  (x)  is  maiorized  bv 
-  m  ' 


Const .m  and  Const. 


for  then  we  have 


-32- 


HR  ( • )  II  ,  <  /  Const  .m*dx  +  / 

(a. 4)  1^1-^ 


j  Const  •  I  ■' I  dx  ^ 


2  3^  7T 

_<  —  Const  .m  +  2 ‘Const  .log  |  x  |  |  ^  ^  Const .  log  m. 


Since  0  <  R(k.)  <1  we  have 


|R  (x)|  ^  2  •  y  (R(k)|  <  4m; 
|k|<2m 


furthermore,  summation  by  parts  yields 


sln(-y)*R  (x) 
2  m 


k]  <2m 


R(k)»[sln(k  +  “  sin(k  -  •i)x]  |  < 

z  ^  — 


<4+  y  |R(k+l)  -  R(k)| *1 sin(k  +  y)x| , 

l<lkT<2m-l 


and  since  R(k)  are  assumed  to  decrease  monotonically 


|r  (x)(  K  -  K  Const.  -T — r  , 

“  |sln(f)|  ~  ^ 


which  completes  the  proof. 


-n- 


B.  THF:  nFCAY  RATE  OF  THE  FOERIER  COEFF  LC  lENTS-RE  V  I  S  ITED 


In  Section  3,  we  concluded  that  the  quantities 
E^(t)  -  II  •  ,  t )  II  2  >  satisfy  for  k  >  m,  the  recursive  inequality 


(3.14) 


L“(x) 


(b.  1) 


2 

Eo,.(t)  <  - max  E  (t)  +  e  ^E  (^). 

2e/k.k-  t/  ,  - 

^T<t 


'2k 


In  this  section  we  complete  the  details  for  the  solution  of  these  recurrence 
relations,  and  obtain  that  for  k  >  2®*m  we  have 


,,  8®E„*N“  s  ,-s+l  ,  2 

(b.2)  E2^(t)<[_^)  .E^  +  (1  +  _JL^  >  -"k  t 


8®E^*n2  s+1 


c/Tc*  k" 


i.e.,  (3.15)  holds.  For  s  =  0,  (b.2)  is  reduced  to  (3.13);  now  assume  that 
(b.2)  is  valid  for  any  k  >  2®*m:  in  particular,  for  k  >  2®''’^*m  we  can 
use  (b.2)  with  k  replaced  by  7  ,  y  >  2®.m,  and  obtain  that 


( b .  3 )  max  Ej^(t  )  <  [ 


s+1 


oS+l„  .,2 

^  8  E„N  s  , -s  ,  2^ 

— •!!„  Ml  +  -  ■_  )  •e-''  nE„. 

e/k»  k  e/k»k ' 


'O' 


Furthermore,  we  have 


8®E-N^  s+1  s  -s  7 

(b.4)  E2^(-^)  <  -E^  +  (1  +  ^  ^  --2-^ 


e/l<«  k 


e/1c«k‘' 


r-) 


'0' 


Fsing  (b.3)  and  (b.4)  to  upper  bound  the  right  hand  side  of  (b.l)  we  find 


-34- 


R  .N 

i  •  ( 


2e/k«k.  e/k*k 


EqN 


2e/k*  k^ 


s+1 

8  E^*N“  s+1 

)  -^0^ 


o  s+ 1 „  2  „ 

8  E-N  s  ,  -s  ,  2 
r,  ^  0  ^  -4  ek  t  „ 

(.1  +  - 5-J  *6  -E 

e/lc*  k 


0 


„  .2^  S^E-N^  s+1 
,  -2ek  t  f  Os 
+  e  .[___]  .E 

e/k*k 


0 


^  ,  2  ft  E„N  s  ^  -s  ,  2 

,  “2ek  t  ^,  .  0  -2*4  ek  t  ^ 

+  e  *[1  +  — z — 

e/k*k 


O' 


s+2 

The  first  of  the  four  terms  on  the  right  is  less  than  ( - — — ^ )  *^0’ 


the  sum  of  the  remaining  three  terms  does  not  exceed 


e/1<*k 


s+1  ,-s,2^ 


1  + 


'0 


e/k*k 


g+  \ 

and  hence  for  k  >  2  'm  we  have 


-4  ek  t 


En. 


E,k(t)  <  ( 


s+2 


e/1?*  k 


T 


)  •£„+(!+ 


s+1  2 

8^  E-*N  s+1  ,-s  ,2^ 

0  -4  ek  t 


•E, 


e/T?*k‘^ 


which  completes  the  induction  proof  of  (b.4), 
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